Coulomb drag at zero temperature 
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We show that the Coulomb drag effect exhibits saturation at small temperatures, when calculated 
to the third order in the interlayer interactions. The zero-temperature transresistance is of the or- 
der h/(e 2 g 3 ), where g is the dimensionless sheet conductance. The effect is therefore the strongest 
in low mobility samples. This behavior should be contrasted with the conventional (second or- 
der) prediction that the transresistance scales as a certain power of temperature and is (almost) 
mobility-independent. The result demonstrates that the zero-temperature drag is not an unambigu- 
ous signature of a strongly-coupled state in double-layer systems. 

PACS numbers: 73.23.-b, 73.50.-h, 73.61.-r 



Coulomb drag effect has proven to be a sensitive probe 
of electron-electron (e-e) interactions. The phenomenon 
is usually observed [H, Q, 0, 0, 0, @| in double-layer sys- 
tems, where electrons interact through the long-ranged 
Coulomb forces. A current, passing through one of the 
layers (the active layer), induces a voltage across the sec- 
ond (passive) layer. The ratio between the two, the so 
called drag transresistance pn, carries a valuable infor- 
mation about the state of electrons in each of the layers, 
as well as the nature of their mutual interactions. 

The transresistance for weakly interacting electrons 
was calculated 0, HI, 0, [ToL 11 | in the second order in 
the screened interlayer interaction and found to be given 
by 



p D {T) = 0.12. 
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where d is the separation between the layers. Ep and 
fei? are the Fermi energy and momentum correspondingly 
and k is the inverse Thomas-Fcrmi screening radius. This 
result is in a reasonable agreement with a number of ex- 
periments P, 0, 0, H[ . Its main feature is the quadratic 
temperature dependence, which may be traced back to 
the phase volume accessible for the interlayer e-e scat- 
tering. The second order effect requires the electron-hole 
asymmetry, i.e. the difference in velocity between elec- 
trons and holes on the opposite sides of the Fermi surface. 
Such an asymmetry scales as Ep 1 for each of the two lay- 
ers, giving rise to the factor Ep 2 in Eq. (JTJ) . The latter 
serves as the dimensional scale, which normalizes the T 2 
dependence. 

On the other hand, the systems with strong interlayer 
correlations are predicted to exhibit a nonzero drag tran- 
sresistance (oc h/e 2 ) even at zero temperature. The exam- 
ples include ID charge density waves at exact commensu- 
rability [l^ |. as well as Quantum Hall bilayer structures 
at the total filling factor v = 1 (TjJ. In the latter system 



the effect was likely observed experimentally in Ref . [14 1 . 
This raises a question if pp> (0) may serve as an unambigu- 
ous indicator of a strongly-correlated state in a system 
at hand. I.e. whether the drag transresistance undergoes 



a quantum phase transition between the weakly-coupled 
state, where it is strictly zero, and a strongly- coupled 
phase, where it is finite. 

In this Letter we give a strong argument against such 
a scenario. We show that pd(0) ^ already in weakly 
interacting bilayer systems. To this end we evaluate the 
transresistance in the third order in the (screened) in- 
terlayer interactions and find a constant temperature- 
independent contribution 



p D (T) = 0.27 
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T < h/r . 
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where g = 25.8fcf2/pn is the dimensionless conductance 
(here pa is the resistance of the single layer) and r is the 
elastic scattering time. Drag effect, saturating at small 
temperatures, is therefore not an automatic indicator of 
a strongly correlated state. 

There are general reasons to expect that the third or- 
der effect may be qualitatively different from the second 
order one, Eq. {!]). Indeed, the third order transresis- 
tance does not rely on the electron-hole asymmetry. This 
is because the corresponding linear-response diagrams 
involve four-leg vertices (see below) which do not van- 
ish within linearized dispersion relation approximation 
(i.e. in the electron- hole symmetric case). Therefore, 
the result is expected to be independent on the Fermi 
energy, Ep. Since we are interested in the lowest tem- 
peratures, it is natural to focus on the diffusive regime, 
where T <C h/r. In this regime there are no any other 
relevant energy, which may provide a scale for a temper- 
ature dependence. Hence, the temperature-independent 
result, Eq. (2j, is not entirely unexpected. Moreover, 
the four-leg vertex, mentioned above, is known to play a 
central role in the low-temperature transport of diffusive 
metals. It is exactly this object that gives rise to singu- 
lar Altshuler-Aronov (AA) corrections to the miralayer 



conductance [15 1 



Coming from another perspective, it is certainly un- 
usual to find a temperature-independent result for the 
quantity which relies on the e-e scattering rate. Indeed, 
the latter is proportional to the available phase volume 



2 



around the Fermi surface, which scales as T 2 . However, 
in addition to the occupation numbers the scattering rate 
is proportional to a certain matrix element (the overlap 
integral of six wave functions for the third order process, 
considered here). In diffusive systems such matrix ele- 
ments are known to be singularly enhanced in the limit 
where all involved states are close in energy [l6j|. It is 
exactly this enhancement that leads to singular e-e in- 
teraction effects in the low-temperature diffusive limit 
fvfy . In the case of the third order transconductance in 
2D the smallness of the phase volume is exactly com- 
pensated by the divergence of the corresponding matrix 
elements. This yields the temperature independent tran- 
sresistance, Eq. (J2j> - The diffusive enhancement of the 
matrix elements is less pronounced in cleaner systems. 
Hence, in the clean limit, g — > oo, the zero temperature 
drag, Eq. ([2]), disappears. 

There are two limitations on the applicability of Eq. ([2|) 
at the very low temperatures, (i) Once the tempera- 
ture length Lt = y/hD/T, where D is the diffusion 
constant, reaches the sample size L, the growth of the 
matrix elements is saturated. As a result, prj oc T 2 at 
T < £rhi here Erh = hD/L 2 is the Thouless energy, 
(ii) If the sample size is very big, one may enter the 
regime of disorder and/or interaction induced localiza- 
tion. The relevant temperature scale is that where AA 
correction a\j = e 2 /h[g — tt~ 1 ln(/i/TV)] [l5| is significant, 
i.e. T ~ (/j/V)e _7r9 . At smaller temperatures the diffu- 
sive approximation breaks down and our result, Eq. ([2]), 
is not applicable. 

A natural question is why in experiments of e.g. 
Refs. [3,1113 

the low-temperature saturation of Pd(T) 
was not observed. In order to answer, one may estimate 
the saturation temperature T* by equating Eqs. {JJ and 
©. This way, one finds T* « E F (k F d)g~ 3 / 2 . Employ- 
ing the parameters of e.g. Ref. Q: Ep « 60K, g w 100 
and fcf d as 4, one finds T* as 0.25K and the residual re- 
sistance pd « 0.4mfi as it follows from the Eq.©. At 
the same time the lowest temperature reported in Ref. [H 
T as 0.5K and the corresponding drag pr> ~ 0.65mf2 are 
just above the expected saturation. The similar situa- 
tion is true regarding most of the other reports of the 
Coulomb drag Refs. 0, BUS]. 

It is rather likely, though, that the saturation observed 
by Lilly et. al Ref. [l8| in v = 1 Quantum Hall bilayer 
system in the composite fcrmion regime is a manifesta- 
tion of Eq. (|2|). Indeed, it was shown [3] that the dif- 
fusive corrections in the composite fcrmion regime arc 
rather similar to those in zero magnetic field. Virtu- 
ally the only difference is a substantial downward renor- 
malization of the composite fermion conductance g c f , as 
compared to the zero field one, g. Estimating g c f as 10 
and E C J as 5K for the samples of Ref. (HI, one finds 
T* as 0.15K and pd(0) « 2fi in a good agreement with 
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FIG. 1: The four- leg vertex, central to the third-order drag 
effect, as well as to the intralayer A A correction. External 
wavy lines represent fluctuating vector potentials; the ladder 
is the diffusion propagator T>(r — r',oj). 



tures are needed. 

The four-leg vertex, which is a building block for di- 
agrams of the third order drag transconductance is de- 
picted in Fig. [JJ It describes an induced non-linear in- 
teraction of electromagnetic fields through excitation of 
electron- hole pairs in a given layer. The vertex is non- 
local because of the diffusive propagation of the electron- 
hole excitations within the layer. The latter is encoded 
in the propagator 
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V a (q,uj) = — 



1 



v a D a q 2 



(3) 



Ref. 18(. To verify Eq. © more experiments in zero 



magnetic field with smaller g or/and smaller tempcra- 



where v a is the density of states of the layer a = 1, 2 and 
D a is its diffusion coefficient. Notice that the dimension- 
less conductance is expressed as g a = v a D a . 

We work with the Keldysh technique 2^, 21 , 2^ ■ In its 
framework the fluctuating electromagnetic potentials ac- 
quire an additional index: classical (cl) or quantum (q), 
which stay for symmetric and antisymmetric combina- 
tions of the fields propagating forward and backward in 
time, correspondingly. The proper indices are indicated 
in Fig. [TJ The fact that the four-leg vertex of this very 
structure is unique in the leading order in l/g a may be 
rigorously proven within Keldysh non-linear sigma-modcl 
|2l| . In fact, it is exactly this vertex which gives rise to 
the singular AA correction [15| . The latter is obtained 
by pairing one classical and one quantum electromagnetic 
potentials, while the two remaining ones represent an ex- 
ternal (classical) electric field along with induced current 

It is convenient to work in a gauge, where the Coulomb 
interactions are mediated by the longitudinal vector po- 
tentials, rather than the scalar potentials. An advantage 
of using such a gauge is that both internal and external 
potentials, as well as the current sources are all expressed 
through the same type of field. This makes the structure 
of the vertex, Fig. [TJ particularly symmetric. Moreover, 
the gauge may be chosen in a way that the propagator 
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of the longitudinal vector potentials V a p = 2i(A^A%) 
automatically includes the vertex renormalization by the 
disorder 



21] 



(D a q 2 - iLu)(Dpq 2 - iuS) ' 



(4) 



where Vj^(q, ui) is the 2x2 matrix of retarded intra and 
interlayer interactions calculated within random phase 
approximation (RPA). The latter is the solution of the 
standard Dyson equation 
where 



MUM V R = V + V flV R , 



V 
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(5) 

Note that the polarization operator n(q,oj) has no off- 
diagonal elements, reflecting the absence of tunneling be- 
tween the layers. 

We are now on the position to evaluate the third order 
drag transconductance. The corresponding diagrams are 
constructed from the two vertices of Fig. [TJ one for each 
of the layers. Fig. @. Remarkably, there are only two 
ways to connect them, using the propagators ((31) (recall 
that (A%A q p) = 0, ;22J). The analytic expression for the 
sum of the two diagrams of Fig. [2] is given by 



a D = 32e 2 Tg 2 g 2 J ^1^2^ Im [z>i (q, lo)V 2 (q, u) 

(6) 
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Via(«, uj)V 12 (|-Q,|-Q) V12 (| + Q, I + f 



The two functions !Fi(iO, 0) and J-%{oj, ^) originate from 
the integration over the fast electronic energy s, Fig. [TJ 
in the active and passive layers correspondingly. In the 
dc limit they are given by 

Fi(w, fi) = T A [fi(n + u/2) - 6(fi - w/2)] , (7a) 



.F 2 (cj, fi) = 2 - 6(0 + w/2) - B(n - w/2) + B(lu), (7b) 



B(aO = £c«Lh(£) 



(7c) 



To make the farther calculations more compact, we re- 
strict ourselves to the identical layers. We shall first con- 
sider the experimentally most relevant case of the long- 
ranged coupling, where d > 1. Here k = 2ire 2 h' is the 
Thomas-Fermi inverse screening radius. In this limit the 
effective interlayer interaction potential, Eqs. (J4j) , ([5j) , ac- 
quires a simple form 



V 12 (q,oj) 
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FIG. 2: Two diagrams for the drag transconductance an in 
the third order in the interlayer interactions, Vi2(g,w), de- 
noted by wavy lines. The intralayer diffusion propagators 
T> a (q,uj), Eq. ([3]), are denoted by ladders. 



Next, we substitute Eqs. ©, J7J) and © into Eq. © and 
perform the energy and momentum integrations. The in- 
spection of the integrals shows that both energies to and 
n are of the order of the temperature u> ~ Q ~ T (in com- 
pliance with the phase volume considerations) [23j. On 
the other hand, the characteristic value of the transferred 
momenta is q ~ Q ~ ^T/(Dnd) <C ^/T/D, cf. Eq. ©. 
Therefore one may disregard Dq 2 as compared to iu> in 
the expressions for T> a (q,Lo), Eq. ([3]), approximating the 
product T>iT>2 in Eq. ^ by — cj~ 2 . This factor represents 
the diffusive enhancement of the matrix elements, men- 
tioned in the introduction. Such spatial scales separation 
implies that the four-leg vertices, Fig. [TJ are effectively 
spatially local, while the three interlayer interaction lines 
are long-ranged. 

Rcscaling energies by T and momenta by y/T/(Dnd), 
one may reduce the expression ([6]) for the transconduc- 
tance to <j£> = (e 2 /h) g~ 1 (K,d)~ 2 x [dimensionless inte- 
gral]. The latter integral does not contain any parame- 
ters and is free from divergences in all directions. It is 
thus simply a number that may be evaluated numerically 
[24] ]. In the limit ajj <C {e 2 /h)g a the transresistance is 
related to ov by po — o Dh 2 / '(e 4 51172) , resulting finally 
in Eq. ©. 

To emphasize the fact that the scale separation, dis- 
cussed above, is not crucial for having the low tem- 
perature saturation, we briefly consider the case of the 
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short-ranged interlayer interactions, V-^(q, u>) = Vq. The 
latter may be relevant, if interactions are screened by 
e.g. metallic back gate. One employs then Eqs. ©, dU) 
and rescales the energies by the temperature, while the 
momenta by y/T/D. This way the transconductance, 
Eq. ^ , once again reduces to the dimensionless and pa- 
rameter free integral. The latter is convergent in all di- 
rections and may be readily evaluated, resulting in 



PD = 0.01 4 4 (^o) 3 



(9) 



Notice, that the effect is expected to have the negative 
sign for the short-ranged attractive interactions. This ob- 
servation may have relevance for oppositely doped double 
layer structures. 

The low temperature saturation of the Coulomb drag 
was discussed previously in Refs. [25| and [26|. Both 
of them considered essentially different and somewhat 
more exotic mechanisms. The zero temperature satura- 
tion suggested in Ref. 25[ relies on the assumption that 
the electrons in both layers are scattered by exactly the 
same disorder potential. Ref. [26| focuses on the strongly 
coupled regime, where the pairing order parameter is sup- 
pressed by disorder. 

To conclude, we studied the Coulomb drag phe- 
nomenon in weakly interacting bilayer systems. We 
found that effect saturates at small temperatures, when 
calculated to the third order in the interlayer interac- 
tions. The saturation of drag relies on the presence of 
disorder and scales inversely with mobility. It does not 
require, though, any correlations of the disorder poten- 
tial in the two layers. The effect was possibly observed 
in Ref. (l8j , although more experiments in lower mobility 
samples and zero magnetic field are highly desirable. 
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where Ko(rp) is the modified Bessel function, which is 
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Here p = \J —\x and fi± = y/— i(±y + x/2), the radius 
is normalized by y/Dnd/(2T) and x = ui/T; y = fi/T. 
Finally the number of interest is given by 
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